Tricritical point in explosive percolation 
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The suitable interpolation between classical percolation and a special variant of explosive percola- 
tion enables the explicit realization of a tricritical percolation point. With high-precision simulations 
of the order parameter and the second moment of the cluster size distribution a fully consistent tri- 
critical scaling scenario emerges yielding the tricritical crossover exponent 1/ipt = 1-8 ± 0.1. 



PACS numbers: 64.60.ah, 89.75.Da, 64.60.al 

An explicit realization of tricritical percolation has 
been a longstanding challenge. The diluted Q-state Potts 
model yields, at the tricritical point, Ising exponents, for 
Q— but a direct construction of this tricritical perco- 
lation has not yet been achieved [Tj. The recent appear- 
ance of explosive percolation, that exhibits a discontinu- 
ous transition, provides the essential new element which 
enables us to find the way of establishing tricritical scal- 
ing as introduced in this Letter [5H1]. 

A discontinuous percolation transition is observed 
when the growth of the largest cluster is systematically 
suppressed [1], promoting the formation of several large 
components that eventually merge in an explosive way 
[5]. Several aggregation models, based on percolation, 
have been developed to achieve this change in the nature 
of the transition [21 lU [BHS] • These models are generally 
classified as explosive percolation, the name given in the 
original work that triggered the field In that work, a 
best-of-two product rule is proposed to occupy the bond 
which minimizes the product of the mass of the merg- 
ing clusters O IT0Hl3] . More recently, this procedure has 
been generalized to a best-of-m product rule in random 
graphs [H] and regular lattices to study its perco- 
lation and transport properties. For different values of 
m, previously reported models are recovered pi [51 [TB]. 
Potential applications are, for instance, the growth dy- 
namics of the human protein homology network [17] and 
the identification of communities in real systems |16| . 

The larger the set of bonds m considered in the prod- 
uct rule, the lower the probability that the occupied bond 
is related to the largest cluster [TS] , promoting the com- 
pactness of the percolation cluster. Inspired by the Q- 
state Potts model which, in the limit Q— ^1 its magnetic 
transition is related with the geometrical transition of 
percolation [18] . we also propose a hybrid model, where 
an additional parameter is included which allows us to 
interpolate between this discontinuous explosive regime 
and the one from classical percolation, characterized by 



a continuous transition. We obtain a nonequilibrium tri- 
critical percolation where explosive percolation is diluted 
with classical percolation. Without dilution the process 
is explosive (discontinuous transition) while under maxi- 
mum dilution classical percolation is recovered (continu- 
ous transition). A multicritical behavior appears due to 
the addition of this new degree of freedom (degree of di- 
lution) and a new set of critical exponents is found at the 
tricritical point [1^1 HO] • We also analyze the size depen- 
dence of the order parameter and the susceptibility when 
a best-of-ten product rule is considered, and numerically 
corroborate the hypothesis of a discontinuous transition. 
In the best-of-m product rule, at each iteration, m bonds 
are uniformly and randomly selected from the set of un- 
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FIG. 1. Size dependence, for the best-of-ten rule {q — 0), of 
the maximum in the susceptibility x and the order parameter 
Poo (inset), at the percolation transition, for square lattices 
with size , with L ranging from 32 to 4096 sites. For large 
system sizes, the order parameter converges to a constant 
value and the susceptibility scales with L'*, where d is the 
dimension of the system. Results have been averaged over 
10'* samples. Error bars are smaller than 0.6%. 



2 



occupied ones. From this set, only the one minimizing 
the product of the mass of the clusters it connects is 
occupied. In the limit of m = 1 the model recovers clas- 
sical percolation [3T], while m = 2 corresponds to the 
Achlioptas process [5]. Nagler et al. [55] have reported 
that for Erdos-Renyi graphs (mean-field limit) a strong 
discontinuous transition solely occurs when m is equal to 
the total number of bonds in the system, while for lower 
values of m the transition is rather weakly discontinuous 
- the gap shrinks with the system size and it is zero in 
the thermodynamic limit. Also Smoluchowski equations 
for another model - similar to the product rule - reveal 
a continuous transition in the mean- field limit |T4] . On a 
square lattice, Monte-Carlo simulations reveal a discon- 
tinuous transition, characterized by some properties of 
a continuous one [31 [T31 US] . Increasing the number of 
considered bonds m promotes the formation of compact 
clusters, delays the percolation threshold, and, above an 
intermediate value, improves the conductivity of the sys- 
tem \TE\ . Here we take the best-of-ten case characterized 
by a percolation transition for the fraction of occupied 
bonds Pc = 0.55975 ± 0.00002, obtained from the cross- 
ing of the wrap probability for different system sizes |24j . 
In Fig. [l] we show the dependence on the system size, 
for this model, of the maximum in the susceptibility (%) 
and the order parameter (Poo) at the percolation thresh- 
old. As in Ref. [1], we define the susceptibility as the 
sum over the square of the finite cluster (i.e., nonperco- 
lation) mass - also known as mean cluster size - and the 
order parameter as the fraction of sites belonging to the 
largest cluster. The susceptibility scales linearly with the 
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FIG. 2. (color online) Order parameter (Poo) as a function 
of the fraction of occupied bonds (p), for different values of 
the degree of dilution q. Two different regimes are found. 
For lower values of q the transition is discontinuous while 
for larger values it is continuous. Increasing q the percolation 
threshold goes from the one for the best-of-ten product rule to 
the one for classical percolation. Results have been averaged 
over 10^ samples for square lattices with lateral size of 4096, 
and the error bars are below 0.2%. 



number of sites in the lattice and the order parameter, 
for larger system sizes, converges to a nonzero value, as 
in Ref. [13] for the best-of-two rule, consistent with the 
hypothesis of a discontinuous transition. Let us now de- 
fine a hybrid model with a new parameter q, representing 
the degree of dilution. With probability q the next oc- 
cupied bond is randomly chosen among the unoccupied 
ones (to = 1), while with probability 1 — q the bond is 
selected following the best-of-ten product rule (to = 10). 
Adding an additional parameter to the model, allows us 
to interpolate between different regimes. In the limit of 
q = 1 we recover classical bond percolation. On the other 
limit, for q = 0, the model boils down to the best-of-ten 
product rule. In Fig. [2] we see the order parameter for 
different values of q. One can clearly identify two differ- 
ent regimes. While for large values of q the transition is 
continuous, for smaller values, the transition is discon- 
tinuous. The larger the value of q the lower the fraction 
of occupied bonds at the percolation threshold. 

The hybrid model is characterized by two control pa- 
rameters, the fraction of occupied bonds p and the degree 
of dilution q. The former is the control parameter that 
triggers the transition between a nonpercolative and a 
percolative system. The latter interpolates between the 
two different regimes, namely, the best-of-ten and the 
classical. Figure [3] shows the two-parameter diagram of 
the model. The horizontal axis is the probability q and 
in the vertical axis the fraction of occupied bonds p. The 
transition lines correspond to the values of the fraction 
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FIG. 3. (color online) Phase diagram of the model. The 
horizontal axis is the degree of dilution q and the vertical 
one the fraction of occupied bonds p. The (green) dashed 
line is the first-order transition line and the (red) solid line 
is the continuous transition. Both lines meet at a tricritical 
point {qt,Pt)- Below the transition lines, there is no infinite 
cluster in the system. Above the transition lines an infinite 
cluster exists. The (blue) arrows stand for the relevant scaling 
fields fiq and fip (discussed in the text). The values have been 
obtained from the crossing of the wrap probability, with error 
bars below 0.2%. 
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of occupied bonds at which the percolation transition oc- 
curs for different values of q. The discontinuous line is 
dashed while the critical line related to a continuous tran- 
sition is solid. All points on the critical line belong to the 
classical percolation universality class. Both lines meet 
at the tricritical point {qt,Pt)- To analyze the crossover 
between the two different regimes we follow the theory 
of Riedel and Wegner ^9\. Two relevant scaling fields 
need to be defined. One tangent to the critical line at 
the tricritical point (/i^) and another perpendicular to it 
(/ip), as depicted in Fig.|3] In this coordinate system the 
critical line is described by [TH US] 



(1) 



where (fit is the tricritical crossover exponent. To 
determine the tricritical point we plot, in Fig. [4j 
the relation between both relevant scaling fields, 
along the critical line, for different values of q 
(0.49,0.50,0.51,0.52 and 0.53). The best straight crit- 
ical line, in a log-log plot, as described by Eq. ([T]), is 
obtained for qt = 0.51 ± 0.01, where the index t stands 
for tricritical. For this value, the inverse of the tricritical 
crossover exponent is l/tpt — 1.8 ±0.1. The scaling fields 
can be obtained from the transformation 
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with 6 ~ 0.022, the angle of fip with the horizontal axis. 
For the typical percolation transition obtained by cross- 
ing the critical line in the direction of the Hp scaling field, 
the tricritical point is characterized by a new set of tri- 
critical exponents - not only the crossover exponent ipt, 
but also, for example, /3t, 74, and Vf To obtain these 
exponents we consider, in Fig. [5j the size dependence of 
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FIG. 4. (color online) Scaling of the critical line, according to 
Eq. ([1]), for different values of q. A straight line is obtained for 
the tricritical value qt — 0.51 ± 0.01. The obtained crossover 
exponent (pt is such that l/ift = 1.8 ± 0.1. 
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FIG. 5. Size dependence of the maximum in the susceptibility 
(x) and for the order parameter (-P00), at the tricritical point 
of the hybrid model. System sizes range from 32 to 4096 
lateral sites. Results have been averaged over 10^ samples. 
Error bars are smaller than the point size (below 0.5%). 



both the maximum in the susceptibility (x) and the order 
parameter (Poo); a-t the tricritical point. The former gives 
= 1.84 ± 0.05 and the latter Pt/i't = 0.09 ± 0.01. 
The scaling ansatz for the susceptibility is given by 



x{p-VuL)^L'^^'-^F\{p~pt)L^'''^ 
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where F [x\ is a scaling function which, for large x scales 
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FIG. 6. (color online) (a) Finite-size scaling for the sus- 
ceptibility (x) for q — qt, according to the scaling ansatz of 
Eq. ID Results for square lattices with lateral size of 1024, 
2048, and 4096 sites, are averaged over 10* samples. The 
value of 7t = 2.9 ± 0.1 can be obtained from the slope of the 
dotted line, (b) Crossover scaling for the susceptibility (x) 
for different values of g = {0.515,0.52,0.53,0.54,0.55} and 
the largest system size (4096^ sites). Dashed and dotted lines 
correspond to the two different regimes (details in the text). 
Results have been averaged over 10'' samples of a square lat- 
tice with 4096 sites of lateral size. 
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as x~^* . Figure |6|a) shows the scahng for three different 
system sizes, namely, 1024, 2048, and 4096, giving 74 = 
2.9 ± 0.1. Also, from Figs. [5] and |6ja), we find ut = 
1.56 ± 0.09 and, consequently, /3f = 0.14 ± 0.03. 

For the complete tricritical crossover scaling of the sus- 
ceptibility, X is considered as a homogeneous function on 
the relevant scaling fields. 



(4) 



where I is a scaling parameter, h is the ghost field of 
percolation, and 6t can be obtained from 5t — {-jt+Pt)/ Pt 
[26] . In the absence of ghost field {h = 0) and with 
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liq ^^"^^ one has. 
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where -Fcross [v] is a scaling function such that, -Fcross [y] ^ 
for y < 1 and F^^ss [y] - y'^' for y » 1 [^D]. Fig- 
ure |6]^b) shows the tricritical crossover scaling for differ- 
ent values of q, close to the tricritical point. Two different 
regimes are clearly identified. For low values of ix^jp^'J'^^ 
(below 1), corresponding to the larger values of g, re- 
sults are consistent with 7 = 43/18 as expected for the 
classical percolation universality class (dashed line). For 
large values, the exponent of the scaling function corre- 
sponds to the one from the tricritical point reported in 
this work. This scaling gives the tricritical crossover be- 
tween the discontinuous (explosive) percolation and the 
continuous one. 

Our model can be further generalized to different val- 
ues of m (number of preselected bonds). As discussed 
previously, the larger the value of m the stronger the 
signs of a discontinuous transition for the pure best-of-m 
model [ini . As a consequence, the level of dilution at the 
tricritical point qt must increase with m. 

In short, we have introduced a diluted tricritical per- 
colation model, where the degree of dilution allows us 
to interpolate between the classical percolation regime, 
characterized by a continuous transition, and the best-of- 
ten one where, on a square lattice, a discontinuous tran- 
sition is observed. We have presented a two-parameter 
diagram for the model to systematize the different per- 
colation regimes. Two transition lines were identified: a 
discontinuous and a critical line. Both lines meet at a 
tricritical point characterized by a new set of tricritical 
exponents different from the ones of the classical per- 
colation universality class. Finally, we also analyze the 
tricritical crossover scaling for the susceptibility. 
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